home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / SSPDI.z / SSPDI
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. SSSSSSSSPPPPDDDDIIII((((3333FFFF))))                                                            SSSSSSSSPPPPDDDDIIII((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      SSPDI   - SSPDI computes the determinant, inertia and inverse of a real
  10.      symmetric matrix using the factors from SSPFA, where the matrix is stored
  11.      in packed form.
  12.  
  13.  
  14. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  15.       SUBROUTINE SSPDI(AP,N,KPVT,DET,INERT,WORK,JOB)
  16.  
  17. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  18.      On Entry
  19.  
  20.      AAAAPPPP REAL (N*(N+1)/2)
  21.         the output from SSPFA.
  22.  
  23.      NNNN INTEGER
  24.         the order of the matrix A.
  25.  
  26.      KKKKPPPPVVVVTTTT INTEGER(N)
  27.         the pivot vector from SSPFA.
  28.  
  29.      WWWWOOOORRRRKKKK REAL(N)
  30.         work vector.  Contents ignored.
  31.  
  32.      JJJJOOOOBBBB INTEGER
  33.         JOB has the decimal expansion  ABC  where
  34.         If  C .NE. 0, the inverse is computed,
  35.         If  B .NE. 0, the determinant is computed,
  36.         If  A .NE. 0, the inertia is computed.
  37.         For example, JOB = 111  gives all three.  On Return Variables not
  38.      requested by JOB are not used.
  39.  
  40.      AAAAPPPP contains the upper triangle of the inverse of
  41.         the original matrix, stored in packed form.
  42.         The columns of the upper triangle are stored
  43.         sequentially in a one-dimensional array.
  44.  
  45.      DDDDEEEETTTT REAL(2)
  46.         determinant of original matrix.
  47.         Determinant = DET(1) * 10.0**DET(2)
  48.         with 1.0 .LE. ABS(DET(1)) .LT. 10.0
  49.         or DET(1) = 0.0.
  50.  
  51.      IIIINNNNEEEERRRRTTTT INTEGER(3)
  52.         the inertia of the original matrix.
  53.         INERT(1)  =  number of positive eigenvalues.
  54.         INERT(2)  =  number of negative eigenvalues.
  55.         INERT(3)  =  number of zero eigenvalues.  Error Condition
  56.  
  57.      AAAA division by zero will occur if the inverse is requested and  SSPCO  has
  58.      set RCOND .EQ. 0.0 or  SSPFA  has set  INFO .NE. 0 .  LINPACK.  This
  59.      version dated 08/14/78 .  James Bunch, Univ. Calif. San Diego, Argonne
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. SSSSSSSSPPPPDDDDIIII((((3333FFFF))))                                                            SSSSSSSSPPPPDDDDIIII((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      Nat. Lab.  Subroutines and Functions BLAS SAXPY,SCOPY,SDOT,SSWAP Fortran
  75.      ABS,IABS,MOD
  76.  
  77.  
  78.  
  79.  
  80.  
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.